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Exact diagonalizations indicate that the effective 1- 
dimensional behavior (sliding Luttinger liquid phase) of 
the frustrated spin-1/2 crossed chain model, predicted by 
Starykh, Singh and Levine [Phys. Rev. Lett. 88, 167203 
(2002)], persists for a wide range of transverse couplings. The 
extension of the other phases (plaquette valence bond crystal 
and Neel long range order) is precised. No clear indication of 
a coexistence of these two phases is found, at varianc e with a 
suggestion of Sachdev and Park ( :ond-mat /01082 14 
PACS numbers: 75.10.Jm; 75.50.Ee; 75.40.-s 



I. INTRODUCTION 

Among the family of frustrated two dimensional spin- 
1/2 antiferromagnetic (AF) models, a crossed chain 
model (CCM), which may also be viewed as an exten- 
sion of the checkerboard antiferromagnet, has recently at- 
tracted interest in two distinct ranges of parameter 
The CCM Hamiltonian reads: 



H — J\ Si.Sj 



(i.i) 



<l,3> 



«i,j>> 



where the exchange J\ couple nearest-neighbor pairs 
< i,j > of spins on a square lattice and exchange J 2 
next-nearest-neighbor pairs << i,j » on a checker- 
board pattern of plaquettes as shown in Fig. [l| Alter- 
natively, J 2 may be viewed as the intra-chain exchange 
between spins on the diagonal chains and J\ as the cou- 
pling constant between these crossed chains. Both J\ and 
J 2 are assumed positive, describing AF couplings. This 
model interpolates between the AF Heisenberg model on 
the square lattice for J2 = and decoupled AF Heisen- 
berg chains for J% = 0. When Ji = J 2 , one has the 
checkerboard antiferromagnet which is a 2d analog of the 
pyrochlorc antiferromagnet. 

The physics of the spin-1/2 model is well established 
at three points. For J2 = it has collinear Neel long 
range order (LRO) with gapless AS Z = 1 elementary ex- 
citations (magnons). At J2 = Ji it is a valence bond 
crystal (VBC) with LRO in singlet plaquettes its el- 
ementary excitations are gapped and carry integer spins 
When J2/J1 = 00, one has the physics of decoupled 
F Heisenberg chains, their elementary excitations are 
gapless spin-1/2 deconfined spinons. 




FIG. 1. Left: the crossed chain magnet. The spins sit at 
the vertices shown by bullets, the coupling constant between 
nearest-neighbor pairs is Ji, the coupling constant on the di- 
agonal bonds of the checkerboard of plaquettes is J2 . ui , 112 
are the unit vectors of the Bravais lattice. Right: the Bril- 
louin zone associated to this lattice. F, ki (i = 2,5), M are 
the six different wave vectors present in the N=36 sample. 



Recently, motivated by experimental results on 
Cs 2 CuCl A § and Na 2 Ti 2 Sb 2 §], Starykh, Singh 
and Levine (SSL) have investigated the CCM in the 
limit J 2 lJ\ S> 1 where it describes "weakly" coupled 
chains Q. They concluded that, the chains behave then 
as quasi-decoupled, realizing a SU(2) "sliding Luttinger 
phase" E^B. Usually the Id AF Heisenberg behavior 
is unstable to non frustrating 2d couplings giving rise to 
Neel LRO or VBC LRO. SSL argue that in the CCM the 
Id behavior should survive for finite non negligible J\/ J 2 
due to the frustrating nature of the transverse coupling. 
Since the pioneering works of Majumdar-Gosh and Hal- 
dane, many ladder models have been studied giving birth 
to various Id behaviors [p| [llf . The case of the spatially 
anisotropic couplings of the triangular lattice typical of 
Cs 2 CuCli, appears in some aspects rather exceptional 
with non negligible transverse couplings and neverthe- 
less an effective Id behavior with a continuum of gapless 
excitations looking very much like Fadeev spinons Jl2| . 

The first purpose of this paper is to investigate SSL 
prediction using exact diagonalizations (ED) and try to 
figure out the possible extension of the quasi-ld regime 
with decreasing J 2 /J\. The second purpose will be to 
study the range of stability of the VBC phase around 
J 2 /Ji = 1. We shall also address a new suggestion 
of Sachdev and Park ||] on the possible coexistence of 
phases with two different order parameters between the 
pure VBC and Neel phases. 
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II. RESULTS 



The ED calculations have been carried out on N = 
16, 32, 36 samples which display all the symmetries of 
the infinite lattice, using periodic boundary conditions 
as described in p| . Additional calculations using twisted 
boundary conditions were also performed on the TV = 16 
sample ( specifically for J2/J1 > 1) which indicated that 
the lowest ground-state energies are obtained for periodic 
boundary conditions: it might be inferred from these re- 
sults that incommcnsuratephases are unlikely in the ab- 
sence of a magnetic field [Q . 

The general evolutions, as a function of J2 {J\ = 1), of 
the ground-state energies per spin E/N and the spin gaps 
A are displayed in Figs, g, || and [|. As shown in Fig. ^ 
the point of "maximum frustration" where E/N has its 
maximum, is reached for J2/ J\ slightly smaller than one, 
close to the checkerboard point. On might notice that 
the finite size evolutions between N=32 and N=36 are 
hardly visible on the scale of these figures for J2/J1 < 1 
but are dramatic for J2/J1 > 1- We will show below that 
this last behavior is related to the quasi- 1-dimensional 
behavior of the model in this range of parameter. 
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FIG. 3. Spin gap A vs J2 (Ji = 1), same symbols as in 
Fig. [| Extra lines show the spin gap values for independent 
chains of length L = 4 (full line) and L — 6 (dashed line) for 
J2 > 1.25. 
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FIG. 4. Ground-state energies per spin E/N vs J2 for 
3% > 1, same symbols as in Fig. H. Extra lines show the E/N 
values for independent chains of length L = 4 (full line) and 
L = 6 (dashed line) for J2 > 1.25. 



FIG. 2. Ground-state energies per spin E/N vs J2 
(Ji = 1) for J 2 < 1.25: N = 16 (tripodes), N = 32 (crosses) 
and N = 36 (full squares). The dotted lines are guides for 
the eye. 



A. Quasi Id behavior 

ED results show a large extension of Id behavior from 
J2I J\ = 00 down to J2/J1 = 1.5 and probably below. 

First, Fig. [| shows that E/N is quite close to the 
energy values of decoupled chains of same length for 
J2/J1 > 1.10. Comparison of results for N=16 and N=32 
shows that increasing the number of chains pushes the 
energy per spin toward the value of uncoupled chains 
(Note that the N=16 sample has 2x2 chains of length 
4, the N=32 sample has 2x4 chains of length 4, whereas 
the N=36 sample has 2x3 chains of length 6). Point- 
ing to the same conclusion, the spin gap (Fig. ||) shows 



dramatic finite size effects between N=32 and N=36 for 
J2/J1 > 1.10. 

Second, Id behavior down to at least J2/J1 = 1.5 is 
supported by spin-spin correlations: 

s(ij) =< Si.Sj > (2.1) 

and dimer-dimer correlations patterns: 

D(ij;kl) =< Sj.SjSfc.S/ > - < S».Sj >< S k .Si > . 

(2.2) 

Values of s(ij) computed in the ground-state of the N — 
36 sample at J2/J1 = 1.5 are displayed in Table |. 

They are largest and decrease slowly on a line of di- 
agonal bonds and are quite small if sites i and j belong 
to different chains, even if i and j are first neighbors. 
Spins on parallel chains are less correlated than those 
on orthogonal chains which may be understood in a per- 
turbative approach: coupling between orthogonal chains 
involves one " Ji" term instead of two for parallel chains. 
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j 


s(ij) 


j 


s(ij) 


36 


-0.4551363204 


4 


-0.0228887132 


15 


0.1970987646 


2 


-0.0156898505 


22 


-0.2193602830 


9 


0.0102951567 


3 


0.0017845218 


16 


0.0029547330 


10 


0.0006638989 


17 


-0.0040816720 


12 


-0.0009191508 







TABLE I. Spin-spin correlations s(ij) —< Sj.Sj > be- 
tween site i = 1 and site j in the exact ground-state of the 
N = 36 sample for J2/J1 = 1.5. The sites are numbered as in 
Fig. p. Other correlations can be deduced using d v symme- 
tries. Correlations along a J2 chain are displayed in the upper 
left corner of the table , those between orthogonal chains are 
in the upper right corner, those between parallel chains in the 
lower part. 
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OM) 


D 


(k,l) 


D 


29 22 


0.0897246082 


34 27 


-0.0015374269 


22 15 


-0.0523861622 


10 3 


-0.0015374269 


31 6 


0.0110558121 


34 29 


-0.0011366550 


35 4 


0.0074567919 


27 22 


0.0005263317 


28 23 


0.0066321724 


15 10 


0.0005263317 


16 11 


0.0066321724 


11 4 


-0.0003686058 


21 16 


-0.0063063310 


35 6 


0.0001626070 


35 30 


-0.0057705738 


16 9 


0.0001126614 


11 6 


-0.0057705738 


28 21 


0.0001126614 


9 4 


-0.0046514788 


35 28 


-0.0000762647 


17 10 


0.0017481736 


33 4 


-0.0000456425 


34 5 


0.0016564513 


34 3 


0.0000268261 


12 5 


-0.0015704399 


10 5 


0.0000036528 



TABLE II. Dimer-dimer correlations D(ij; kl) (Eq. \j.2j ) 
between the pair of sites (i, j) = (1, 36) and pairs of sites (k, I) 
in the exact ground-state of the N = 36 sample (J2 = 1.5, 
Ji = 1). The values of D for pairs (k,l) on diagonal bonds 
coupled with J2 (left column and upper right corner above 
the line and Fig. a) are all larger than on non-diagonal bonds. 
The sites are numbered as in Fig. H. 




FIG. 5. Dimer-dimer correlations D(ij ; kl) in the exact 
ground-state of the N = 36 sample (J2/J1 = 1.5) between 
reference pair (1,36) and pair (k,l). Full (D > 0) or dashed 
(D < 0)) lines have width proportional to the square root of 
\D\ (a choice done to ease visualization). 



The behavior of the dimer-dimer correlations is iden- 
tical ( Table |ll] and Fig. ||). D(ij; kl) values between the 
pair on a diagonal bond and pairs (k, I) are also 

largest if and (k,l) are on the same chain. They 

decrease by an order of magnitude on nearby regions of 
crossed chains and by two orders of magnitude on parallel 
chains. 

Third, the evolution of the spin-gap A with increas- 
ing size, is not inconsistent with a vanishing value for 
N — ► 00, at least down to J%j J\ = 1.5. In Fig. we dis- 
play the values of A/J 2 vs 1/L for N = 16,32 (L = 4), 
N = 36 (L = 6) and compare them with the spin-gaps 
of the AF Heisenberg chain. At J%j J\ — 5 the values of 
A/J 2 (not shown) would be indistinguishable from the 
Heisenberg chain values and only slightly deviate from 
the latter down to J2/J1 = 1.5. The extrapolation for 




FIG. 6. Ratios A/J2 of the spin-gaps to the intra-chain 
coupling J2 of the chains at J2/J1 = 2 (open losanges), 
J2/J1 = 1.5 (open squares) and J2/J1 = 1-25 (open trian- 
gles). The spin-gaps of the AF Heisenberg chain are shown 
by crosses. All N — 16 values are very close to the L — 4 value 
of the Heisenberg chain. The dotted line joins the N = 32 
and N = 36 results for J2/J1 = 1.5. 



L — ► 00 at J2/J1 = 1-5 is problematic due to system 
sizes: nevertheless a linear extrapolation of the two sets 
of data for N — 32 and TV = 36 points to a zero gap: look- 
ing to the sign of the deviation from the 1 /L asymptotic 
limit (which is negative) one is indulged to conclude that 
the system is still gapless in the thermodynamic limit for 
J2/J1 = 1-5. This could even hold at J2/J1 — 1.25, but 
the extrapolation is then more uncertain with present 
sizes. 

Fourth, a lengthy but straightforward calculation 
shows that the quantum numbers of the ground-state and 
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first excitations of the CCM model for J%fJ\ = 5 (Ta- 
ble III) are directly related to the quantum numbers of 
the ground-state and first excitations of chains of length 
L = 4 or 6 (respectively for the N — 32 and 36 sam- 
ples). The effect of the transverse inter-chain coupling 
may be considered in a first approach as a minor quan- 
titative effect. The lowest S=l excitations consist of a 
set of 4 (N=16), 8 (N=32), respectively 6 (N=36) quasi 
degenerate states. Their parentage to the first excitation 
(labelled a in Table III) of the unique chain can be ex- 
actly traced back. Their wave vectors are vectors of the 
side of the square Brillouin zone (points k$ and M of the 
Brillouin zone in Fig. [I] for N = 36) . Appropriate su- 
perpositions of these quasi degenerate states create S=l 
excitations which are essentially localized on one given 
chain! 




FIG. 8. Same as Fig. § at J 2 /Ji = 1.5. 




FIG. 7. Spectrum of eigen-energies E in each total sector 
of spin S for TV = 36 at J2/J1 — 5. This spectrum is quite 
close to the uncoupled chains one (see Table ^ ). The lowest 
energies E(S) in each S sector increase linearly with S up to 
S — 3 as shown by the dashed line. 

The quasi degenerate states with total energy ~ 77—78 
in Fig. are associated to the creation of excitations 
with 2 "a-quanta" on different chains (levels 2a in Ta- 
ble |iTj| ) . This process gives birth to states with total spin 
S = 0, 1, 2. Endly we should notice that the linear in- 
crease with S of the lowest eigen-states in each spin sector 
is a last proof that one can add to the system up to 3 ex- 
citations of spin 1 (Fig. |^) which are essentially identical 
and have extremely weak interactions. 

Decreasing the ratio J2/J1, the spectra testify of 
stronger interactions between multi-spinons excitations. 
In particular, a new feature does appear: a splitting of 
the singlet 2a excitations with a downward shift in en- 
ergy of some of them which reach the level of the la- 
excitations at J2/J1 = 1.5 (Table |ffl|, Fig. |and Fig. |). 
One can see in Fig. ^|b this set of singlet states with 
excitation energy AE ~ O.6J2 well separated from the 
other singlets. This is a new feature if we compare it to 
Fig. ^c for J2/J1 = 5, where the spectrum is merely that 
of uncoupled chains. But it is by no means related to 
the physics of the neighboring VBC phase as can be seen 
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FIG. 9. Lowest excitation energies above the ground-state 
AE/ J2 vs wave vector k for J 2 = 1 (a) J2 = 1.5 (b) J2 = 5 
(c) with spin S = (crosses) , S = 1 (tirets) , S = 2 (open 
squares) . 



by comparing this spectrum with the VBC spectrum on 
the same sample Fig. ^a. As explained above the la- 
excitations can be described as purely " unidimensional" , 
but these 2a-levels are excitations which cannot be lo- 
calized on a single chain. The large energy gained with 
inter-chain coupling suggest some interrelation between 
chains. Yet the spin-spin and dimer-dimer correlations 
in these excitations are still strongly anisotropic as could 
be seen in Tables |fv|, [v| and Fig. [h]. Note that the 2a- 
excitations occur at wave vectors which can be obtained 
by combining wave vectors of the la-excitations on the 
sides of the Brillouin zone. In the thermodynamic limit 
the 2a excitations are expected to cover densely the Bril- 
louin zone. We conjecture that the appearance of these 
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0.591242101 
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1.246710497 
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0.635033794 





k 3 




1 












2a 


1.246757039 
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0.640143484 
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1.246757966 
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2a 


0.645292084 
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b 


1.302775637 








-1 


b 


1.266511222 





M 


-1 


-1 


b 


0.839974642 
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-1 


-1 
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1.288135982 





fe 5 




-1 


2a 


1.000230994 
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k 3 




-1 
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1.309791780 
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1 


-1 


2a 


1.001481120 
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1.309841802 
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k 3 




-1 


2a 


1.008659631 
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-1 


1 












2a 


1.309846926 
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k,4 




1 


2a 


1.029661098 
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-1 


-1 












2a 


1.309863513 
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-1 


1 


2a 


1.032523222 
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-1 












2a 


1.349565139 
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-1 


-1 


2a 


1.037137527 
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1 


-1 


c 


1.521999231 


1 


2tt/3 




2a 


1.349566020 





k.4 




-1 


c 


1.060994796 


1 


A-2 




-1 


(1 


1.802775637 


1 


tt/3 




2a 


1.349574718 
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1 


-1 
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1.100373750 
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-1 
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1.352732556 
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1.119357730 
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1.352736474 
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-1 


2a 


1.120391023 
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A~4 




-1 












2a 


1.359019058 
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r 


-1 


-1 


c 


1.135860677 


1 


A-2 




1 












2a 


1.359020101 
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k,4 




-1 


c 


1.142923268 


1 


A-2 




-1 












2a 


1.309863513 
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1 


-1 


3a 


1.143097244 
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A;r> 




1 












2a 


1.399283498 


2 






1 


3a 


1.143396461 


1 


M 


i 






2.802775637 


2 


7T 


1 


2a 


1.399292379 


2 


A'4 




1 


c 


1.145902142 


1 


At, 




-1 




2.802775637 





tt/3 




2a 


1.399293748 


2 


r 


1 


1 


2a 


1.155611046 


1 


k 3 




-1 



TABLE III. Energy of excitations of the simple Heisenberg chain of length 6, of the CCM for Ji = 1 and Ji = 5 (respectively 
1.5) on the N = 36 sample. For each eigen-state we display the total spin S, the k eigen-vector (as described in Fig. [l]), the 
characters of the wave function respectively in a n/2 rotation (R) and in a mirror symmetry (a). The index n indicates how 
the excitations of the 2d model can be build from excitations of the Id Heisenberg chain. 



excitations could be the precursor of a collective mode 
supposed to exist in a crossed sliding Luttinger liquid 




j 


B(ij) 


j 


s(ij) 


36 


-0.4044253289 


4 


-0.0200188490 


15 


0.1666623287 


2 


-0.0403867716 


22 


-0.1505704767 


9 


0.0164521956 


3 


0.0505168767 


16 


-0.0129298755 


10 


-0.0536890650 


17 


0.0523125280 


12 


-0.0470765522 







TABLE IV. Spin-spin correlations in the 4th excited state 
of the spectrum at J2/J1 = 1.5. This state is one of the 
excited states of the 2-dimensional mode. As described in 
Table | it is a singlet state of the trivial RI (k = (0, 0), S = 0, 
R= 1, <y = 1). 



All these results (spectra, correlations, scaling of the 
spin gap) indicate a Id behavior which extend down to 
J2/J1 — 1-5, i.e. well beyond the weak coupling limit 
investigated by (SSL). The 2-dimensional S—0 excitation 
mode singled out for Jij J\ = 1.5 on the N=36 sample 
may be an illustration of the collective mode of the " 



OM) 


D 


(k,l) 


D 


29 22 


0.0421174442 


35 30 


-0.0043780206 


22 15 


-0.0099685219 


11 6 


-0.0043780208 


31 6 


0.0345888875 


9 4 


-0.0070999014 


35 4 


0.0055568217 


17 10 


0.0017481736 


28 23 


-0.0011288331 


34 5 


-0.0004575768 


16 11 


-0.0011288331 


12 5 


-0.0037470380 


21 16 


0.0139907353 







TABLE V. Dimer-dimer correlations in the 4th excited 
state of the N = 36 sample (J2/J1 = 1.5). 

crossed sliding Luttinger liquid" : it is highly anisotropic 
as could be expected. Its detailed properties should nev- 
ertheless be taken with some care due to the numerous 
constraints associated with periodic boundary conditions 
on these small samples. 

The precise location of the lower limit of quasi Id be- 
havior is clearly out of reach of exact ED calculations. 
Nevertheless major modifications of the spectra when 
J2/J1 varies are a good qualitative indication of this tran- 
sition. 

At the isotropic point Jij J\ = 1 the ground-state is a 
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FIG. 10. Same as Fig. |5| in the 4th excited state of the 
N=36 sample (J 2 /Ji = 1.5). 



limit. This is a set of low lying states appearing in each 
total spin sector 5 up to 5 ~ "v~N, belonging to specific 
IRs of the space group, which collapse into a degenerate 
ground-state as 1/N for N — > oo and evolve as 5(5 + 1) 
with 5. For collincar Neel LRO there is one QDJS for 
every 5. Presently these states belong to the trivial IR of 
the space group for 5 even and to the k = (0, 0) IR, even 
under a mirror symmetry and odd under a tt/2 rotation 
representation, when S is odd. In a plot of the spectrum 
vs 5(5+1) the QDJS form a line well separated from the 
macroscopic excitations (magnons) which only collapse 
as 1/VN as shown in Fig. for J 2 = 0. The ~ 1/N 
collapse of the QDJS implies that the spin-gap A also 
closes as ~ 1/N. 

On the other hand the beginning of the plaquette phase 
may be estimated by the appearance of the quasi degen- 
erate S=0 doublet (Svbc) sketched in the previous sub- 
section. 



VBC with LRO in 4-spin S=0 plaquettes f4|: this phase 
does not break SU(2) symmetry (it has a spin gap), but 
breaks the space symmetry of the lattice. This last fea- 
ture manifests itself by the existence of a very low ly- 
ing 5 = excited state which collapses to the ground- 
state in the thermodynamic limit. The absolute ground- 
state belong to the trivial irreducible representation of 
the space group (wave vector k = (0,0), even in all op- 
erations of Ci V ) and the first excited S=0 state, which 
embodies the VBC symmetry breaking is a state with 
wave- vector (tt,tv), odd under a 7r/2 rotation around O 
(These two states form the doublet called Svbc m the fol- 
lowing). When J2/J1 increases beyond 1, the first man- 
ifestation of the instability of the VBC phase appears 
in the modification of the symmetries of these low lying 
singlets, which is concomitant with the appearance of 
quasi 1-dimensional excitations in the S=l sector: these 
two manifestations appear as soon as J2/J1 = 1.05 (resp 
1.10) for N = 32 (resp 36). The transition shift to larger 
J2/J1 values with increasing sizes. Its location in the 
thermodynamic limit is difficult with present sizes. At 
J2/J1 = 1.5, as seen above, all characteristics (spectra 
and correlations functions) point to a sliding Luttinger 
liquid. At J2/J1 — 1.25, the scaling of the spin-gap (and 
the spectra) also suggest this behavior. The transition 
might occurs just below this point. 



B. Coexistence of Neel LRO and VBC order 

We now turn to the investigation of the extension of 
the Neel and the VBC phases and the question of their 
possible coexistence. 

As discussed in JT^JT^] the simplest way to monitor the 
extension of Neel LRO is to search for the so-called quasi 
degenerated joint states (QDJS) which allow to break 
SU(2) and lattice space symmetry in the thermodynamic 




S(S+1) 

FIG. 11. Spectrum at J2 = (Heisenberg model on the 
square lattice) vs S(S + 1) for N = 36. The QDJS character- 
istic of collinear Neel LRO at the bottom of each S are well 
aligned (dashed line) and clearly separated from the other 
lowest excitations (see text). 

The destruction of the QDJS and the appearance of 
Svbc are not necessarily simultaneous. Both the Neel 
and the VBC phase may in principle co-exist in a gapless 
phase. In a recent paper Sachdev and Park , following 
earlier investigations in Ref jl7j] , have suggested that the 
destabilization of collinear Neel LRO by 2th order tran- 
sition into a VBC phase must generically occur within 
such a coexistence region. 

Linear spin-wave calculations HH predict that Neel 
LRO survive quantum fluctuations till J2/J1 ~ 0.75. Up 
to this value the characteristic features of Neel LRO are 
rather well observed in the spectra and in the evolution of 
A and E/N with N. Well below J 2 /J\ = 0.75, however, 
one may notice a slight curvature of the line of QDJS and 
there is small irregularity in the evolution of A and E/N 
with TV: A is slightly larger and E/N slightly lower for 
N = 36 than for N = 32 as might be seen in Fig. || and 
Fig. ||. Nevertheless A seems to extrapolate rather well 
to zero at N — ► 00 at J2/J1 = 0.65. And for J2/J1 = 
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FIG. 12. Spin gaps A (full squares), separation of the 
two singlet states of £vbc (open squares), gap between the 
ground-state and the third singlet (crosses) at J2/J1 = 0.85. 



0.75 the 1/N extrapolation of the N — 16,32,36 values 
of A point to a spin-gap at most very small, ~ 0.06, 
a value of the order of the corrections to the asymptotic 
regime which are expected to be large and neg ative as we 
approach the boundary of the Neel phase |p^|l9| . So the 
Neel phase might persist up to this point. On the other 
hand, up to Jij J\ = 0.75, the separation of the two £vbc 
states appears likely to remain finite as N — > 00. Thus 
the VBC phase is not reached. 

By contrast results at J2/J1 = 0.85 point to a finite 
spin gap and a degeneracy of the two states of Svbc 
which remain separated from the other singlet states for 
N -v 00 (Fig. [if). 

These results lead to conclude to a transition between 
the two phases close to J2/J1 ~ 0.75 in agreement with 
the spin-waves result and that a coexistence region of the 
two phases either occurs at most in a very small range or 
is absent. 



III. SUMMARY 

Exact diagonalization have allowed to characterize 
three phases of the CCM: the (tt, tt) Neel phase at small 
diagonal couplings, 4-spin S=0 plaquette VBC at the 
isotropic point, and sliding Luttinger liquid for J2/J1 > 
1.25. The transition between the Neel phase and the 
VBC phase occurs at J2/J1 ~ 0.75. No evidence of a 
coexistence region predicted by Sachdev and Park || is 
found. The CCM exhibits Id behavior for a large range 
of parameters (J2/J1 greater than w 1.25) where the 
transverse coupling between the chains is of the order 
of the intra-chain coupling. This Id behavior is likely 
to be truly realizing the SU(2) "Luttinger sliding phase" 
with deconfined spinons predicted by Starykh, Singh and 
Levine [jj]. 
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